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The theoretical treatment of surface friction of liquids 
or gases at a solid wall encounters serious difficulties as 
soon as the processes are no longer defined hy the viscosity 
of the fluid alone, but also involve the forces of inertia 
with the probable exception of the flow phenomena in capil- 
laries and the problems of lubricant friction — as is the 
case in nearly all practical problems. All the same, two 
substantial advances have been achieved in this domain with- 
in the last decades; namely, by Prandtl's "boundary layer 
theory" and Blasius' confirmation of the previously sus- 
pected nature of the friction loss in smooth pipes. 

Unfortunately, the results of Prandtl's theory have 
remained confined to a comparatively narrow range, first 
for the more obvious reason, that the paper work involved 
for specific cases is enormous, but then also because its 
physical range of validity is, like the theory of pure 
friction flow in pipes, restricted to narrow limits. Just 
as the pure friction flow, the so-called laminar flow in 
pipes, is replaced by a "turbulent flow" at higher velocities, 
so the laminar boundary layer is replaced by a "turbulent" 
boundary layer. 

The present report deals, first with the theory of the 
laminar friction flow, where the basic concepts of Prandtl's 
boundary layer theory are represented from mathematical and 
physical points of view, and a method is indicated by means 
of which even more complicated cases can be treated with 
simple mathematical means, at least approximately. An 
attempt is also made to secure a basis for the computation 
of the turbulent friction by means of formulas through which 


1 ** 

"Uber laminare und turbulent Reibung." Z.f.a.M.M., 
Vol. 1, No. 4, Aug. 1921, pp . 233-252. 
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the empirical laws of the turbulent pipe resistance can he 
applied to other problems on friction drag. 


MATHEMATICAL IMPOST OE THE BOUNDARY LAYER THEORY 1 

The .problem is restricted to two-dimensional flows; the 
axis y f 0 is chosen as fixed boundary to which the fluid 
adheres . 

The differential equations of two-dimensional flow with 
friction can, by introduction of the stream function \j/ by 
means of the formula 


d\Jj _ d\|/ 

U ^ 

oy ox 

(u and v velocity components in x— and y— direction) 
and elimination of the pressure, be expressed by the single 
equat ion 

SA\!z 3\i/ 3A\|; 0\j/ 3 A\j/ 

— — + =« r A Ay ( 1 ) 

9 1 3 y 3x 3 x 8y 

>3 ^ 2 

o o q, 

where A denotes the operation A = — — + - — — , v = — 

dx e oy p 

kinematic viscosity ( q, = viscosity, p = density of fluid). 
The boundary layer theory refers to flow phenomena for which 
at same distance from the wall the friction shall exert no 
perceptible effect on the velocity field, so that for great 
values of y the stream function changes into an assumedly 
known potential function (x, y, t). At the wall itself 

both velocity components u and v are to disappear. In 
order to meet both conditions, first put 



It is clear that for small values of y the first two 
terms cancel out, ieaving only the stream function ± 

( the str eam function of the boundary layer flow). This is 

1 A list of references on boundary layer theory is given 
in reference 1 . 
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then so defined that at the wall 
other hand, when 




= 0. 


x> 


by bx 
represents a small quantity, 


On the 

y 

</p~ 


be- 


comes very great for all values of y differing appreciably 
from zero; hence, to comply with the first condition — tran- 
sition to potential flow — it is sufficient to determine 

a*, / Won for „ = JL . „ . 

V by Jy - 0 


in such a way that Jv -< — - \ \ j r— 

v 9y \ by J y = o jv 

Thus it is apparent that within the boundary layer (p = 
finite) the first two terms, outside the boundary layer ( P 
very great) the last two terms annul each other. (Strictly 
taken, the u component of the boundary layer becomes the 
u — velocity of the potential flow; for the v component 
the boundary layer flow gives a quantity of the order of magnitude 

JT, which is not contained in the potential flow.) 

Introduce formula (?) in equation (l), arrange in powers 


of and retain only "the highest terms with 


the introduction of 
b'^ o 


expans ion of 


dy 


p = ■— 
and of 


as variable instead of 


Th cn 
y , the 


b^jp 

bx 




+ ( 




by v by 'y = 

0 

v dy 8 'y 

= 0 


S\J/ 0 

( 




Cl X 


'bxby' y 

= 0 

and lastly, considering 

that 

= 

: 0, 

af f or ds 




_1_ 

b 3 'ii 1 

- . , - -J- ... - .. . 

d'i'l 


S 3 \J/i 

J* 

iHdp 2 dp dxdp 3 

dx 

dp 3 

dp 3 

and, after one integration: 




b b^i 1 5 2 \|/ ! 

d\J/j 




dtdp dp dxdp 

dx 

dp® 

dp 3 


according to the formulas 


P-/^, 


' b x 


= 0 




= 0, 


y * o 


(?) 


= f ( X , t ) 
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or with = u, 

on 

introduced again. 



v, the variable y = r\Jx> being 


du 
o' t 


du 

u - — + 
d x 


du 

dy 


d s u 

dy E 


= f ( x , t ) 


(4) 


in agreement with the Prandtl equations. 

The function f (x, t) is determined by the condition for 

y = 0. Since u must change to 
obtained 


u. 


/ ^ o\ 

\ dy / 


y = 0 


there is 


du 0 du 0 

— - + Uq 

dt dx 


f ( x , t ) 


(4a) 


for r, = oo . 


The friction less potential flow follows Bernoulli's 
equation differentiated along the boundary as streamline 
(p Q = pressure along the wall) 


du 0 du 0 _ l dp 0 

dt 0 dx p dx 

hence 

du du du 1 dp d u 

— + u — - + v — = — — — + v — — 

dt dx dy p dx dy 2 


(4b) 


( 4c ) 


The significance of (4) and (4b) obviously is that the 
assumedly known pressure distribution p 0 along the wall 

which arises from the potential flow is to a certain extent 
regarded as impressed field (of force) for the boundary layer 
flow; the pressure differences perpendicular to the wall with- 
in the boundary layer being disregarded. It is this very 
essential hypothesis in Prandtl' s theory that leads to the 
reduction of the number of equations and the arrangement of 
the entire problem. 
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THE MOMENTUM THEOREM OF THE BOUNDARY-LAYER THEORY 


To 'brine out the physical sense of the boundary layer 
theory the evidence contained in the foregoing equations is 
formulated as follows: 

(a) A boundary layer thickness 8 ( as a function of x) 
is to exist such that for y ^ 8 there no perceptible de- 
viation occurs in the flow pattern relative to the potential 
flow; especially the x— component u of the velocity can be 
put equal to the wall velocity of the potential flow u o 
for y = 8 ( x ) . 


(b) Within the boundary layer itself the pressure is 
only dependent on x and equal to the pressure that corre- 
sponds to the potential flow along the wall.. 

By virtue of the two assumptions (a) and (b) the momen- 
tum theorem in the x direction can be applied to a fluid 
volume bounded by the wall, a short piece of the line 
y = 6 (x) and two cross sections perpendicular to the wall 
at x ar+d x + dx (fig. 1.) The increase of momentum is 
equated to the resultant of the outside forces, which involve 
the pressure difference, and the friction R at the wall as 
outside forces. Since for y = 6 the flow changes into 
frictionless potential flow, the friction at the transitional 
area between the boundary layer and the outer field can be 
ignored . 


Hence 


6 

bn b 

t%J puiy * t* 

0 0 

5 

& p 

bt 


2 , d 

pu dy — u 0 


a* J 


' b p 

pudy * — 8 — — — R 
ox 


(5) 


pudy is the time rate of change of the momentum con— 

6 

b r 2 

talned in the considered volume; -r— / pu dy is the excess 

o x 


Ji 


of momentum leaving the front surface over the amount of 

6 

b n 

momentum entering at the rear; -g— / pudy is the inflow 



SAGA. TM No. 109 2 


6 


so that 
this vol— 


In laminar flow the frictional force R referred to 
unit surface is R = T 0 = p, ^ ^ ( T = shearing stress 

in the fluid.) Later on, it is shown that equation (5) is 
practical also for turbulent flow conditions, if u and p 
are regarded as the average values of the velocity with 
respect to time and To is expressed by a corresponding 
empirical formula. 

Eauation (5) can, of course, be derived' also by in- 
tegration with respect to y from equation (4) with due 
regard to (4a) and (4b), It obviously yields, on the basis 
of plausible assumptions for the velocity profile u (y) 
in the boundary layer (0 < y < 5), simply a differential 
equation for 8, that is, for the boundary layer thickness 
as a function of x and t. Limited to stationary processes 
it affords an ordinary differential eauation of the first 
order for 6 as a function of x, so that the development 
of the boundary layer can be followed by comparatively 
simple calculations. The subsequent report by K. Pohlhausen 
(reference l) contains the calculations for a number of 
practically important cases, so that this method need not 
be gone into further. His calculations show that in all 
cases computed by Prandtl’s partial differential equations 
the approximate method ensures results commensurate for all 
practical purposes. In this manner a further development 
of the theory is made possible even where the solution of 
the partial differential equations is extremely tedious, if 
not impos s ibl e , 


volume per unit length of the side area y = 8 (x), 

£ J> 

u — — / pudv indicates the momentum entrained with 

4 

ume of fluid. 


LAM 111 AH AND TURBULENT BOUNDARY LAYER 


The simplest and practically most important case that 
the boundary layer theory deals with is the frictional re- 
sistance of a plate towed in a fluid at rest parallel to 
its own plane. Taking the case of two-dimensional motion 
and referring the motion to the assumedly static plate, the 
problem is as follows: The parallel flow with the uniform 
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velocity U is given as a potential motion, the frictional 
boundary is to start in the origin of the coordinates z = 
y = 0 and for x ^ 0 he given by the axis y = 0. The 
boundary layer thickness and the wall friction is to be 
computed as a function of x. This problem has already been 
solved by Blasius (reference 2); he found that the boundary 

layer thickness increases with J~k . Computing the friction 
drag for a plate of length l and width 1 yields the fric- 
tional force (for friction on both sides) 


W 


= 1.327 ,/npiir 


( 6 ) 


or, if put, as usual, 


\'i = c f y F 


U 

2g 


(7) 


TJ 

where the resistance is referred to the velocity head — , 

2g 

the surface F and the specific weight of the fluid ‘Y = pg, ! 


W = 1.3 27 


/ — TT 

/JL T1 L 

Ul 2g 


( 8 ) 


The coefficient of the frictional drag c^ is a function 
of the Reynolds number, or "reduced velocity" R, in case the 
nondimensional auantity: velocity x plate length divided by 
coefficient of kinematic viscosity is introduced as such, so 
t hat 


c_ = 1.327 — (8a) 

f yir 

Blasius indicated, in a later report (reference 3) based 
upon measurements, that formulas (8) and (8a) are no longer 
valid for large Reynolds numbers, that rather a sudden change 
occurs in the nature of the resistance and presumably in the 
state of flow, similar to that occurring in pipe flow at the 
critical limit. On the other side of the sudden change the 
resistance increases at more than the 3/2 power of the velocity 
hence the resistance coefficient in equation (7) decreases 
s 1 ower than — — . 
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Next, it is assumed that the laminar boundary layer, for 
which the Pr andt 1— Bias ius theory gives the afore-mentioned 
results, is replaced by a "turbulent boundary layer," in 
which — as for the turbulent flow in pipes — the velocity is 
subjected to continuous fluctuations in magnitude and direc- 
tion. The first consectuence of the fluctuations — when 
plotting the streamlines of the average flow — is that the 
shearing stress is not caused by the sliding of the adjacent 
fluid portions alone; the portion of the shearing stress 
corresponding to the friction becomes small relative to the 
momentum transport owing to the irregular convection of the 
supplementary velocities. Up to now, it has not succeeded 
to explore in some way the nature of this momentum convection 
apparently obeying statistical laws — and to make the fluc- 
tuation phenomena accompanying the turbulent flow amenable 
to a theoretical study. In this respect the present article 
contributes nothing to the solution of the puzzle. The task 
undertaken here merely involves the introduction of plausible 
assumptions for the distribution of the average values of the 
velocity within the boundary layer, which are based on the 
empirical law of turbulent motion in pipes and the applica- 
tion of the previously derived momentum equation to the 
equilibrium of the boundary layer. It results in relations 
of the turbulent friction at a towed plate which are in very 
good agreement with experience. 


THE TURBULENT PLOW IN SMOOTH PIPES 


The laws of flow resistance in pipes have been the sub- 
ject of an unusually large number of experiments. But the 
empirical material has not improved much up to within recent 
date because the different degrees of wall roughness had been 
frequently ignored and the tests were not referred to the 
physically correct parameter, that is, the Reynolds number, 
(The only resistance formula, so far, which allows for the 
relative roughness and the Reynolds number is that by R. von 
Mises (reference 4.) In many instances no consideration was 
given to the fact that the constant velocity profile in the 
pipe is formed only after a fairly long "convection path." 
Bias ius merits the credit of having found an empirical 
formula by analyzing the material and comparing the best ex- 
periments for smooth pipes which very accurately reproduces 
the nature of the flow resistance over a, wide range. Accord- 
ing to it the pressure drop for a circular pipe is referred 

y ® 

to the velocity head of the average velocity 
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l u 2 fv 

h * X , X = 0.316 / — (9 ) 

d 2g V vd 

where l is pipe length, and d pipe diameter. 

By this formula, which represents the experiments ex- 
tremely well over a wide velocity range the pressure drop is 
proportional to the 7/4 power of the average velocity as 
against the previously held, conception that the resistance 
lav; above the critical velocity would approach the square 
law fairly soon. R. von Mises incidentally conjectured that 
contemporarily with increasing velocity, the velocity dis- 
tribution over the section becomes consistently more uniform 
so that the measured parabolic velocity profiles merely 
form a transitional phenomenon and that the profile varies 
continuously with increasing velocity (reference 4.) In 
the technical literature a parabolic distribution independ- 
ent of the velocity is for the most part tacitly assumed. 

The writer concurs with Von Mises to the extent of assuming 
a distribution varying with the Reynolds number but with the 
difference of assuming a well defined distribution function 
as asymptotic form rather than the uniform distribution, 
which the velocity distribution approaches at larere Reynolds 
numbers and on perfectly smooth walls. Hence, the assump- 
tion that in the turbulent as in the laminar zone, at least, 
for large Reynolds numbers for which the resistance law (9) 
holds true, a similar remaining velocity distribution over 
the cross section exists, so that for increasing throughflow 
volume all velocities increase in proportion. Prandtl raised 
the question whether conclusions could be drawn from the em- 
pirical law (9) regarding this velocity distribution. He 
found on the basis of a dimensional analysis that, under 
certain plausible assumptions, the resistance law definitely 
defines the distribution of the velocity in the direct 
vicinity of the wall. The suggestion for the following anal- 
ysis goes back to a conversation with Prandtl in the fall 
of 1920. The publication is with his consent, though the 
process of derivation is somewhat different from his. 

Consider a pipe of circular cross section. If the 
velocity in the pipe axis (f = 0) is indicated by u max , 

the assumption of a velocity profile independent of the 
throughflow' volume and increasing similarly implies that the 

ratio -p" — is a definite function of — only (r = distance 
max a 

from pipe axis, a = pipe radius). 
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Hence, the first assumption reads; the velocity at dis- 
tance r from the pipe axis can be put 

u = u max =p(|) 

/V \ 

cp (■§;•) being independent of u Das . On doubling the velocity 
in the center all velocities are doubled. 

The second assumption states; the velocity distribution 
in vicinity of the wall, that is, near r = a, is to depend, 
aside from the physical constants p and p, only on the 
distance from the wall T| = a — r, and further, on the 
shearing stress (frictional force) T n transferred to the 

wall. Hence, for small values of p 

U = f ( u,, p , T q , p ) (ll) 

Specifically, the auantity u is to be independent of 
the pipe dimensions, that is, of a, for small values of 
p. This assumption is based upon the plausible concept 
that the velocity distribution next to the wall is inde- 
pendent of the other boundaries of the flow, so that a 
definite relation exists between the friction on a wall ele- 
ment and the immediately adjacent velocity distribution. 
Visualize eauation (ll) developed by increasing powers of p; 
the first term of the development to read 


u = f x ( u, p , t q )p x 


( 11a) 


x to be defined later. 

The third assumption contains the empirical resistance 
law: on doubling the velocity the pressure drop and the 

shearing stress at the wall T Q is to be increased as 

1 : 2 7 / 4 . 

The dimensional eauality of the left- and right side 
of (lla) can obviously be maintained only when f contains 
the quantities u, p, T 0 also only in powers; for, on 

AT 

bearing in mind that ,/ — and — have the dimensions of 

p n 
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velocities, it is readily apparent that the only possible 
dimensionally correct combination is 

u = b (i°) (-n -Y (12) 

\ p / ^ v ' 

where B is a nond in ens ional constant. 

On the other hand, since u increases according to ( 1 0 ) 
in proportion to the throughflow volume; whereas, T 0 ac- 
cording to the resistance law increases with the 7/4 power 
of the through flow volume, the relation 

i--~ = 4/7, x = 1/7 ( 13 ) 

2 


must apply. 


The first term of a development of the velocity is thus 
obtained as a function of the wall distance 


u = 




( 12a) 


or, with u ( T] ) denoting the velocity distribution in prox- 
imity of the wall, the shearing stress 



( 12b) 


B is a universal constant valid for smooth walls the 
magnitude of which is obviously contingent upon the statisti- 
cal law of the turbulent fluctuation eouilibrium. 


It is somewhat surprising at first to find the differ- 
ential Quotient at the wall to be infinitely great. Since 
no momentum convection can occur on a smooth wall because 
both velocity components disappear, the shearing force must 


be equal to the frictional force 



This expression 


should be infinite according to eolation (12a). The matter 
is explained, however, by the fact that the equations (12) 
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and (12a) must be regarded as an asymptotic expression for 
the velocity distribution at infinitely large Reynolds num- 
bers just as the power law for the flow resistance represents 
an asymptotic law for absolutely smooth walls and for very 
large Reynolds numbers. The true velocity distribution is 
obtained by drawing a tangent with finite slope, say, to 

the velocity curve, so that T = u -g~. (Compare, for ex- 
ample, the interesting measurements in reference 5.) It is 
readily apparent that the point of contact of this tangent 
is shifted to the point n = 0 with increasing Reynolds 
number. But it appears that equation ( 1 2 ) itself represents 
the velocity distribution with sufficient accuracy for moder- 
ate Reynolds numbers. 

The best experiments on the velocity distribution in a 
circular pipe were undoubtedly those made by T. 1. Stanton 
(reference 6), first, because he originated the use of very 
fine pitot tubes in velocity measurements, and second, he 
employed a very long straight entrance section ahead of the 
test section, thus ensuring that the measurements fell in 
the zone where the velocity profile no longer varied per- 
ceptibly. Figure 2 shows Stanton’s velocity values (ratio 
of local velocity to pipe axis velocity) against the wall 
distance, both on a logarithmic scale. It is seen that — 
apart from the first test point, 0.25 millimeter from the 
wall, so that the indication of the pitot tube of 0.33 milli- 
meter in diameter no longer seems reliable — that the test 
points lie very accurately on a straight line of l/7 slope. 1 

For the further applications quantity B in (12a) and 
(12b), which according to the assumptions for smooth sur- 
faces signifies a universal constant of the turbulent flow 
regime, must be determined next. For this purpose it is 
really necessary to know the total velocity distribution from 
the wall vicinity to the center of the pipe, whereas the for- 
mulas (l2a) and (l2b) are valid, for the present, only in wall 
proximity. The chosen method of calculation included the use 
of several appropriate interpolation formulas which satis- 
factorily reproduce the velocity distribution, as measured 
by various experimenters, and change to eauation ( 12a) at the 
wall . 

x It is to be noted that Christen proposed a velocity 
distribution formula according to which the velocity is pro- 
portional to the l/8 rather than l/7 power of the distance 
from the wall (reference 7), A detailed presentation of the 
various distribution formulas is found in Forchheimer ' s work 
(reference 8) as well as in G-umbel's report (reference 9). 
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(a) An extreme case occurs when the formulas which in- 
dicate the velocity is proportional to the l/7 power of the 
distance from the wall are continued to the pipe center. 
Iher ef or e 


/ 

' a - r ' 

\ 1/7 

f y r 

y 

u = u max \ 

^ 3. * 

J ~ u max < 

V 1 a 

) 


( 13 ) 


or for wall proximity 

v p J 


( 1\ 1/7 

V u / 


_1 /7 

n 

= u max ~T7T 
a 


Considering the relation 


do £ . rf 

rra = 2uaT n and T n = Y ——a 


dx 


2 l 


existing "between the pressure drop and the wall stress, the 
calculation of the flow resistance "by Blasius’s formula gives 


Xu ' 
8 


4 /? 


) 


1 

177 


U 


m a x 


a 


1 7 ? 


The ratio of average velocity u, occurring in Blasius’i 
formula, to maximum velocity is, by ( 13 ) , : 


= 0.816 


u. 


max 


Thus with X = 0.316 
B is v vd/ 




1 / 4 


the value for the constant 


l/7 ( — ) 0.816 = 8.57 

V 0.316 / 


0.316 


(b) A better approximation to the measurements is 
afforded in the case where the velocity profile at the pipe 
center is "rounded off" a little. This is best obtained by 
the formula 
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u = 



( 13a) 


wh ere the exponent n can "be chosen arbitrarily; n _ i 
obviously leads back to (13). Figure 3 contains a number 
of measurements of different experimenters along with the 
three curves for n — 1, 1,25, and 2. The test points lie 

almost without exception between curves n = 1 and n = 2, 
Repeating the above calculating process with n = 1.25 and 
n ~ 2, the constant B amounts to 


B = 8.52 with n= 1.25 

B = 8.82 with n= 2.00 

average velocity 

The values of the ratio are 0.838 

maximum velocity 

for n = 1,25 and 0.875 for n = 2. The most reliable 
measurements give 0.84. 1 From this it is concluded that ( 13 ) 
with n ?= 1.25 to 2 represents the conditions fairly accu- 
rately, so that hereinafter B = 8.7 is generally used. 


Thus eouation ( 12a ) must be written 



( 14a) 


If the shearing stress T 0 is expressed as function of 
the velocity, equation (12b) reads 

T o = -774 P lim ) 

' T) as 0 VuTl/ J 

With the values of B obtained on the basis of the 
three interpolation formulas, the enuation would read 


x Also worthy of mention are the measurements by G. J. 
Williams (reference 10) and Gumbel (reference 9), where 
the proportionality factor deer eases a little with increas- 
ing Reynolds numbers and then approaches the limiting value 
0,811. This would favor the simple interpolation formula 
under (a) . But substantially higher values (up to 0.8^) also 
occur, where the effect of entrance length and roughness have 
not yet been fully explained. 



SAC A TM So. 1092 


15 


T 


0.0225 p lim 

ri = 



( 145 ) 


as the general expression for the wall friction in case the 
velocity distribution u (p) is known in the vicinity of 
the wall. The constant in equat ion (l4b) amounts to 0.023? 
for the velocity distribution in the pipe by (14;, as against 
0.0231 and 0.0221 with p = 1.25 and P - 2, by (l4a). 


APPLICATION TO HEAT TRAN SEER 


Jor comparing these formulas with the representations 
expressing the turbulent friction by an apparent increase 
in friction coefficient (reference ll) the shearing stress 
transmitted in a layer distant p from the wall is 


T 


g( P . U , 


P . 


u) 



( 15) 


If u as a function of p and the pressure gradient 
in the pipe are known, the function g can be explicitly 
calculated. Near the wall T must become T 0 . 

du 

Considering (14a) and especially the relation: 7 p = 

u, yields for g: 


g (p, u. p) = 0.8 Op p f ^ 


3/7 


\ p / 


1/7 6/7 

u P 


and putting T = T 0 — as follows from the condition of 


equilibrium for the circular pipe, gives 


g ( P , u . P ) 


0.805 



1/7 e /7 

v Y 


( 15a) 


Y being solely a function of p which becomes p for small 
p. The relation (15a) is applicable to any cross section if 

it is assumed that the ratio of shearing stresses -=f~ is 

1 0 

independent of the velocity and is only a function of the 
lo cat i on . 
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Quantity g is a kind of "turbulent friction coeffi- 
cient" or better expressed a "turbulence factor." 

The portion of true friction at larger Reynolds numbers, 
up to an extremely thin layer at the wall, is vanishingly 
small, hence the shearing stress is to be regarded almost 
exclusively as an average value of the momentum convection. 
This identification is of interest because it makes it 
possible to develop further the analogy between frictional 
resistance and heat transfer in turbulent flow, discovered 
by Reynolds (reference 12) and Prandtl (reference 13), 
Assuming that the momentum transport and the heat transfer 
is accomplished by the same mechanism of the irregular molar 
fluctuating motion, evidently results in two analogous formu- 
las for the shearing force transmitted perpendicular to the 
flow by "turbulent momentum conduction" per unit of surface 
and for the heat volume transferred by "turbulent heat con- 
duction;" 


T = 


0,8 05 



J 6/7 d( 0U) 

v y _-u 

dy 


a =? 0.8 05 



1/7 .,6/7 
v Y 


d( coQ ) 
dy 


( 15a) 


where c = specific heat, 9 * temperature, hence c9 = heat 
content per unit ma.ss. Formula (15a) may be continued up to 
the wall with good approximation if the same pr opor t ional ity , 
assumed for the mechanism of the "turbulent momentum and heat 
transfer," exists for the transfer of molecular momentum and 
heat, that is, for the laminar internal friction and for the 
true heat conduction. As previously pointed out by Prandtl, 
this is evidenced by the fact that for the respective fluid. 


c* 1 1 

the relation — * 1 exists between heat conduction A, 

\ 

friction coefficient u , and specific heat c- This condition is 
approximately complied with, in gases. If differs very 

fS. 

much from unity, as, for example, for water, the formula may 
be extended only to the boundary of the laminar layer next 
to the wall; while the effect of this layer — as will be 
explained elsewhere — can be expressed by a limiting con- 
dition. 
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The formula (15a) enables the heat transfer to be com- 
puted. in all cases, where the "velocity field" of the turbu- 
lent flow is applicable to the average values in time, and 
hence g is known. With this formula, H, Latzko (reference 
14) worked out a number of technically important cases of 
heat transfer to turbulent flows. It succeeded, in particu- 
lar, in showing, that it is incorrect to speak of a "heat 
transfer factor," as is customary in engineering, that the 
heat transfer is rather conditional upon the total arrange- 
ment, It also succeeded in explaining the effect of the 
individual factors and so to organize the occasionally con- 
tradictory experimental material. In this respect the calcu- 
lating possibilities of heat transfer processes appear sub- 
stantially extended beyond the Prandtl analogy conclusions, 
since for the latter a complete agreement in velocity and 
temperature field had to be assumed, while the formula of the 
present report makes the differences between both also amena- 
ble to calculation. 


TURBULENT BOUNDARY LAYER ON THE ELAT PLATE 


The subseouent calculations are based on the previously 
derived equations (14a) and ( 14b ) according to which the 
velocity distribution as a function of the wall distance is 

— Cf) 


•Tf 


i /? 


( 14a) 


if t o, the shearing stress transferred to the wall is 
given, while the shearing stress T 0 is 

p / v \ 1 / 4 

T q = 0.0225 pu ( K ~J (14b) 


u(p) being the velocity distribution in the neighborhood 
of the wall. To apply these relations to the "turbulent 
boundary layer" requires a corresponding formula for the 
velocity distribution. With 5 = boundary layer thickness, 

U = velocity in undisturbed flow, and y = distance from the 
wall, the elementary formula reads 


u 


( 16 ) 
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Equating ( 16 ) to (14a) obviously gives 



that is, the shearing stress 


T 


o 


0.0225 



( 17) 


Equation (17) yields the formula that must be used in 
the momentum equation of the boundary layer as the frictional 
force in order to obtain a theory of the turbulent boundary 
layer which is to replace the Pr andt 1— Bias ius theory for 
the laminar boundary layer. 

Placing, in fact, equation (17) in equation (5) gives 




8 

^pudy = 0.0225 pU E 



Determination of the integrals 


means of (16) gives the differential 
layer thickness 



eouation of the boundary 


_L 1® 

72 dx 


0.0225 



The solution of this equation reads 


4/5 
x 

or, for the length l 


-(t) 


4 / 5 


( 0.0225 ) 


4/5 


1 / 5 




0.37 


l 



1/5 


( 18 ) 


( 18a) 
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The laminar boundary layer grows proportional to y x , 

4/5 

the turbulent layer proportional to x , according to eaua— 
t ion ( 18 ) . 

Now the frictional resistance of a plate of length l 
can be calculated, either by integration of the frictional 
forces along the plate or by applying the momentum theorem 
to the end section for X 35 l . The resistance (both sides) 
follows as 


v x i/5 


fi = — pU 5 1 *» 0.036 pU l (jj-j J 

Referring the resistance through the formula 

U s 


( 19 ) 


W — cfFy 


2g 


to the velocity head, gives the resistance coefficient Cf 


c 


f 


0.072 


R 


0,2 


the Reynolds number R being put at R 


U l 1 
v 


( 19a) 


Figure 4 contains the test data by Gibbons and Wies elsberger 
(reference 15) on comparatively smooth plates, the line for 
c f according to (19a), the Reynolds number and the resist- 
ance coefficient are given on a logarithmic scale. The agree- 
ment is exceptionally good. Gebers (reference 16) obtained 
a slightly higher exponent. It is suspected that towing of 
very long plates is accompanied by inevitable vibrations 
which permit the resistance to in cre a se rap i dly. _ 


^■According to a communication b j letter, Rrandtl pos- 
sessed formula (19a) before the writer did. He indicates 
(cf. Ergebnisse, vol. I) a similar formula with a supple- 
mentary term which allows for the possible existence of 
laminar flow at the front edge of a suitably sharp plate. 
After determination of the numerical factors from older tests 

by Gebers, he gives the formula c^ = 0.073 ^ 1600 


R 


0 • 2 


where the numerical factor in the second term generally de- 
pends upon the degree of sharpening, and should be practically 
vanishingly small for a rounded— off leading edge. 
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An equally good confirmation is afforded by the velocity 
measurements in the vicinity of a towed hoard. Thus the 
points in figure 5 represent the measured velocity distri- 
bution perpendicular to a board towed in water as a function 
of the wall distance and specifically in a section 8.56 meters 
behind the front edee ( reference 1?) . The solid line gives 
the velocity distribution according to equation (16), the 
boundary layer thickness was computed by equation ( 18 ) . 

A comparison of the test data with the curve according to 
equation ( 16 ) discloses* above all, that it is in no way 
necessary to assume a velocity jump at the wall, as com- 
monly reported in the technical literature. The present 
formulas rather represent the rapid decrease in velocity 
next to the wall by the variation of the power curve with 
the exponent l/7 correctly and unrestrictedly. 


LAMINAR FLOW ON A ROTATING- DISK 


As a further illustrative example of applying the meth- 
ods obtained for the calculation of laminar and turbulent 
frictional resistances, the case of a uniformly rotating 
disk is to be analyzed. The laminar state of flow caused 
by a rotating flat disk is of special interest for the rea^- 
son that it represents one of the rare cases in which the 
differential equations of the viscous fluids can be inte- 
grated without ommissions. It offers an immediate check 
on the accuracy with which PrandtJ's boundary layer actua- 
tions yield an approximation. 

The problem is posed as follows: 

The half space x > 0 shall be filled with liauid. 

The boundary plane x = 0 rotates about the x— axis with 
the uniform rotational speed w. What is the state of 
motion in the half space x > 0 with consideration to the 
fluid f r ict ion ? 

Introducing cylindrical coordinates r , • d , x and de- 
noting with c r , c t , c x the velocity components in radial, 

tangential and axial direction, and with p the hydro- 
static pressure, the differential equation of the flow in 
cylindrical coordinates — when, as follows from reasons of 
symmetry, all velocities are independent of d — read: 
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ftr 

ft c 1 
ftr 

ft c . 


+ c 

X 


4- r* ' 
"T 


+ C. 


ftx 

ftct 

ftx 

ft 


i ftp r ft 

+ V 1 

p ftr ^ 


: ♦ -5- (Sr) 

ftr a ftr ^ r ^ 


s 2 Vft 
ft c, 

ftr 8 


r^ct 


1 c t c r .. . 

+ ss 1/ \ + 

(•ftr 8 ftr 


r 

1 ftp 


s ( 21 s ) • s2 °tl 


+ v 


f ft 2 c. 

l; 


r ftr A ftx p ft x ^ ftr ‘ 

and the equation of continuity 


ft c_ c_ Sc 

— L + ~± + x 


1 ftc, 

r ftr 


ftx 8 J 

ftx' 


■J 


= 0 


ftr 


ftx 


\ ( 20 ) 


J 


( 20 a ) 


The construction of the eouations shows that the sys- 
tem (20) and (20a) can he satisfied by the formula 


= rf(x), c j. = rg(x), c = h(x), p = p(x) 


( 21) 


thus yielding for the three functions f, g, h the or- 
dinary simultaneous differential equations 


df 

d 8 f 

. . dg 

d 8 g 

dh 

f - - g* + h — = 

» , 

2fg + h — = 

u , 

* J — 

dx 

dx 8 

dx 

dx 2 

dx 


+ 2f 


( 22 ) 


while the equation 


dh 

dx 


1 dp 

— + y 

p dx 


2 

d h 
dx s 


(27) 


arising from the third equation of the system ( 20) defines 
the pressure distribution p(x). 

Since the fluid is to possess no rotation at infinity, 
but is to adhere to the rotating wall for x ~ 0 , the sys- 
tem of the boundary conditions reads 

f(0) = 0 

g( 0 ) = U) 

h( 0 ) = 0 


f(») = 0 

g( 00 ) = 0 
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The function h( x ) has a finite limiting value for 
x = o°. This means that there is a steady inflow against 
the rotating wall as is to he expected for reasons of 
continuity. Owing to the adherence of the fluid, the ro- 
tating wall acts like a kind of centrifugal fan. Next to 
the wall the fluid is cont inuoxis ly carried to the outside, 
to he replaced hy axial inflow. 

Nor nond im eps ional representation 



( 34) 


is introduced as independent variable and in place of 
f, g, h the functions f, g, h 



id u) 


(34a) 


so that (22) becomes 


f 2 ~ £ S + h 


df d‘f 
d| d | 2 


dg d dh 

?fg + h = r , + 2f = 0 (22a) 

“ di d^ s de 


with the boundary conditions 


f = 0, | = 1, h = 0 for |=0 
f = 0, £ = 0 for |=oo 


so that the equations are independent of all special data 
of the problem. The similarity laws of the problem are 
readily apparent. Since g( £ ) indicates the proportion- 
ality factor of the rotative speed at distance x-ivr 


from the wall to speed of rotation iu , it is clear that 
with increasing velocity only one layer at the wall mani- 
fests perceptible rotational speeds, which decrease with 

increasing velocity and decreasing viscosity as 


On the other hand, it follows from the last of equation 
(24a) that the axial inflow velocity increases at infinity 

as Jim. 
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The equation system (22) can he solved by any numeri- 
cal method or by series expansion. However, it is pre- 
ferred to apply P oh lhaus en ' s method described in reference 
1, for a first approximation. 

It is assumed that the function f and e at a 
distance 6 from the wall are already very little differ- 
ent from zero. From ( 21 ) it follows that the "boundary 
layer thickness" in the present case is constant along 
the wall, hence 6 is independent of | . 

Integration of the first two eouations of (22a) be- 
tween l - 0 and x = A, that is, t = S - £ o , fives 


J> ( 25 ) 


to 


to 


(l £ ~ £ S ) 


df 
h — at 
~ dp 


df 


1 


d£ Jo 


,o |o 

'or 4 g /\ 4g - d £ 

? £i d t i h d | j 


to 

!0 


Partial integration of the second integral, while 

bearing in mind that according to the last eouation of the 
dh 

system (2?a) can be replaced by — 2f , finally affords 


to 

r f 3d t 



df 

II 


to 


4 


f gd| 


m 

L dt Jo 


( 26) 


Put 


1 = a h ( x - C 1 * P fc) 

£ ■ l / s ( 2 + T?) C 1 - jx) 



as approximate expressions for f and g, where a sig- 
nifies a constant that is to be determin"ed T 


It was borne in mind that 
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f = 0, £ = 1 for £ * 0 


df dg 

t ~ ~T> i. ~ ~7 for I = 5 

d£ d| 

and also, as is readily seen from (22a), that 
d £ f 


* - 1. 


d E £ 

,2 

di 


= 0 for 


0 


Numerical calculation of the integrals contained in 
( 26 ) gives 


Co 


f 3 d| = |o [0,0301a 2 - 0,00326a + 0.00159] 




y-’ 

r\ 


g°d£ = to 0,235V 


) ( 28 ) 


t c 


fgd| = io [0,0607a - 0.C0567] 


which, entered in (26), gives two ordinary equations for 
a and |o 


0,09 03a 2 - 0,00978a - 0.23093 = - 
0.2428a - 0.02328 = 


a 

17 ® 

3 

2 |o S 


V ( 29 ) 


J 


The numerical solution gives 


a = 1,026, |o = 2,58 (20) 

On the basis of these data the boundary layer thick- 
ness 5 and the axial inflow velocity c K at infinity 
can then be computed. Obviously 
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6 ^ | 0 y~ = 2.58 y~, - o m = a) 2fd£ = 0.7 08 /iui (3l) 

o' 

Assume, for example, air as fluid with v - 0,14 square 
centimeter per second and an rpm = 600 per minute, that is, 

#5 

uj - _ILS = 62.8 per second: then the boundary layer thickness, 
60 

according to ( 3 1 ) , would he 

8 = 0.122 cubic centimeter 

and the axial inflow velocity 


c^ = — 7.6 centimeters per second 


The most important problem is the calculation of the 
frictional resistance. Assuming the wall bounded by r = a, 
the case is obviously that of a rotating disk with radius a. 
However , the fact that the outer parts of the plane x = Q 
are missing, cannot be without some effect on the motion of 
the fluid, although it is to be presumed that this effect 
remains insignificant, when the thickness of the boundary layer 
relative to disk radius is very small, ap is almost always t he case in practice. 
On these assumptions, the moment of the shearing forces act- 
ing on the disk is simply integrated from r = 0 to r = a, 
or what amounts to the same thing, the angular momentum leav- 
ing in unit time with the fluid at the cylindrical surface 
r = a is computed and equated to the moment of the frictional 
forces. The latter process is preferred. The angular momen- 
tum of the fluid leaving at the cylindrical surface in unit 
time is 

6 


D = 2rra p / c„c_dx = M 


( 32) 


or, by (2l), (24), and (24a) 


n 


M = 2rra 4 pu> 3//2 v 1/<a / fe&i 


J " 


n 


The integral / fgd£ has already been computed. 
o"k 

ing its value from (28) gives 


( 3 2a) 


Enter- 
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M = 0 .9 2a 4 pv 1 ^ 8 w 3/a (33) 

or 

H = 1.84a 3 — (33a) 

2g 

Ua 

where U = au> is the circumferential velocity and R - —rr 

u u 

is the Reynolds number. Equation (33) must, of course, he 
used twice, to obtain the resistance of the disk exposed to 
flow from both sides. 


FRICTIONAL RESISTANCE OF ROTATING DISK IN TURBULENT FLUID MOTION 

The relation governing (33): a frictional moment propor- 

tional to the 3/2 power of the rotational sp?eed at higher rota- 
tional velocities is not borne out in practice. On the con- 
trary, a substantially quicker increase in frictional moment 
is recorded with the rotational velocity. So the assumption 
is made again as for the towed flat plate that a turbulent 
boundary layer is involved, and an attempt is made to secure 
an approximate value for the boundary layer thickness and the 
frictional resistance by applying the momentum equation. 

For the rotating disk two equilibrium conditions are re- 
quired, one in the radial, the other in the tangential direc- 
tion. 


Employing the same notation as in the preceedina section 

and adding T and t (fig. 6) for the frictionsl forces 
r t 

per unit surface at the wall, the momentum Quantities in the 
radial direction are: 


(a) Excess of outgoing momentum Quantity at the cylin- 
drical surface (r + dr)d (for an arc element of opening 
angle l) over to the incoming momentum quantity at area rS 


6 



(b) The radial component of a respectively ingoing and 
outgoing momentum quantity at the front surface (equal to 
the centrifugal force of the rotating fluid volume) 
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o 


These momentum quantities must "be in equilibrium with 
the shearing force T r rdr* hence 


P 


d 

dr 


8 


8 



n 

/ 


/ 


■j 

o 



s 

di 


T r r 


( 34) 


In tangential direction the difference of the turning 
moment of the momentum leaving at the cylindrical surface 
2tt( r + dr ) 8 and entering at the surface 2nrd can be com- 
puted and equated to the turning moment of the frictional 
forces acting on the circular surface. 


2tt 


iL 

dr 


8 



c r c t 4x f 


.1 = 


T t 2rrr2 


o 


(35) 


The formulas for the velocity distribution are according 
to the results of the section Turbulent Flow in Smooth Pipes. 


1/7 


c r = c o u 


-!)■ 


c t = rU) 


(!) 


1 / 7 


(36) 


with due consideration that 


for x = 0 c = 0 , Cj. ~ ru>; for x = 5 c = c. = 0 
r o r t 

Now the integrals in (34) and (35) can be evaluated: 

8 8 


c r dx = 0.207 c 0 6 


c r ctdx = 0.06S1 ruuc 0 8, = 0.02J8 r S oj 2 6 

° ( 37 ) 
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Further, put, in conformity with the assumptions on the 
measure of turbulent friction eouation (14b), 


T = 0.0225 p 


7/1 1 /4 

C U f" 

~T7T 


1 + 


2 -i3 / 6 


V, J 


T t = 0.0225 p(rcu) 


7/4 V> 


1 / 4 


, 1 / 4 


1 + 

/ c 0\ 3 ’ 

__ 

N. ru)/ _ 


3 /e 


by combining the velocity components at the wall and applying 
the friction formula to the resultant. 

With this equations (34) and (35) give the two differential 
equations 


d f 1 

■( o on ' 7 ~ 3 ^ 


dr *- 


207 c s r6 - 0.027S r 3 ui 2 6 = - 0.0225 c 3 r 

0 J ° K Cr,8 


v x 1/4 r 
0 


1 + 


L 



a -.3 /P 


c o y 


d 

dr 


0.06S1 r 3 uuc 0 8 = 0.0225 r^ui v 


! 


4 3/4 b \ 

1,1 ( ) 

v hju64 


1 /4 _ 


/ON 

1 + (-) 

v ru) y 


3/8 


The equations are satisfied if the relationship between 
boundary layer thichness and axial distance r is put as 


c 0 = arm l v 
6 = fSr 3/sJ 


(38) 


and which gives two ordinary equations for a and 3 analo- 
gous to the equation system (29). 


The equations 


0,7456 0,^3 - 0.0278 3 


0.3133 a3 



0 .0225 



( 1 + a s ) 


3 / e 


(? 9 ) 


I 
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given after division 

1,0589 a 2 - 0.0278 = 0 


The numerical solution gives 


p \ 

a ~ 0.162 and hence 8 ~ 0,462 ( — ) 

U) / 


1 / 5 


( 40) 


With these figures the boundary layer thickness becomes 


/ P ^1/5 

6 = 0.462 r ( --- 1 
\r “to 


J 


The section modulus can be computed by*the method given 
in the preceding chapter or else based on equation (55): 


6 

g /"> 5g/P\l/5 

M = 2rra p /' c r c t dx - 0.03S4 a i» p (— =— 

\a w 


y 


( 41) 


and the friction for both sides of the disk 


n / j- 

M = 0,0728 aVp (—\ (41a) 

In accord with the calculations on the towing resistance 
of plates all frictional resistances are then referred to the 
velocity square and velocity head, respectively. With U = 
circumferential speed of the disk, the moment is 


u ^ 

K = 0.146 V — a 3 

2g 


v Ua // 


1/5 


( 42) 


or the resistance coefficient c f as a function of the 

Ua 


Reynolds number of the disk R = 


c« = 0,146 - 

i 5 


p 

1 

yr 


( 43 ) 
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Figure 7 shows Cf plotted, against the Reynolds numbers 

R — both on a logarithmic scale. The plot also contains the 
resistance coefficient 


Cf 


3.68 


( 44) 


secured from the calculation of the laminar boundary layer 
by equation (33a). 

The experimental data were taken from a recently pub- 
lished report on frictional resistance of smooth disks in water 
by W, Schmidt (reference 18). The agreement is good. Of par- 
ticular interest is the fact that the measurements at smaller 
Reynolds numbers fall exactly in the transitional zone be- 
tween laminar to turbulent flow, 1 


NOTES ON ROUGHNESS 


While for perfectly smooth pipes the Blasius resistance 
law is apparently applicable over a wide range of velocities, 
so that it seems more than an interpolation formula, pipes 
with rough sides soon exhibit after exceeding the critical 
point an approximately square relationship between gradient 
and velocity. For this state the pressure gradient may be 
put at 



W 2g d 


(45) 


The experiments by Odell (Engineering, vol. 77, 1904, 
p, 33 and by A, Stodola (steam turbines, 4t.h ed. , Berlin 1910, 
pp , 120—129) on the friction of rotating disks in air give 20 


to 30 percent higher values and a more rapid increase in fric— 

2 + e 

t ion with the circumferential velocity. Odell shows ~ 


with e a small positive digit, Stodola ou 1 * 9 instead of 

u>, * ). Odell’s tests are certainly doubtful, because the 

paper disks which he used, flutter and thus simulate greater 
frictional resistance. In Stodola’s report the higher ex- 
ponent appears to correspond to the roughness of the disk. 
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€ 

X is a function of the relative roughness — ; e denotes a 

quantity with the dimension of a length, which, to a certain 
extent measures the average increase in the wall roughness. 

The ratio of this quantity to the pipe diameter is termed 
the "relative roughness" (reference 4). 

The square law for rough walls is plausible for the 
reason that the frictional resistance is visualized as being 
built up from the individual resistances of the wall pro- 
tuberances, which singly obey the square law. The mechanism 
of frictional resistance is in these cases obviously caused 
by uniform shedding of vortices of well-defined intensity 
and dimensions, as is the case for flows on resistance bodies. 

The flow resistance in perfectly smooth pipes might be 
visualized such that in this instance vortices of dissimilar 
magnitude are separated and float at random ' in the turbulent 
flow, the frequency of the vortices of different intensity 
and size being controlled by some unknown, statistical law. 

By this concept the frictional resistance in smooth 
pipes can be regarded as a fictitious combination of resist- 
ances that correspond to the individual kinds of vortices. 
Assuming that a relationship exists between size of vortices 
and roughness, it may be said that the frictional resistance 
in smooth pipes can be obtained by superposition of the in- 
dividual resistances observed on rough pipes and increasing 
with the square of the velocity, if the individual squared 
resistances are entered with correct weights in the calcula- 
tion. 


It is not without interest that conclusions can be drawn 
about the form of the function X ^ ^ on the basis of this 

conception, and specifically, without knowing the law of 
frequency and the weight function of the individual resist- 
ances . 


In particular, it can be shown that, if the Blasius law 


for smooth pipes holds true, the function 

£\3 / 7 

d 


the formula X 


must have 
e 


at least for small values of — , 

d 


X 0 denoting a constant 
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On superposing the squared resistances by (45) on assump- 
tion of a weight function <p ( e ) , the resistance law for smooth 
p ip es reads 1 : 

00 



( 46) 


Next, it is assumed that the function cp( e ) for small 
values of e compared to pipe diameter is dependent only on 
the physical constants and the velocity distribution in the 
immediate proximity of the wall elements. Specifically, 
cp(e ) is to be independent of‘ the pipe diameter. On the other 
hand, the velocity distribution directly adjacent to the wall 
is according to earlier assumptions entirely contingent upon 
the shearing stress at the particular wall element. There- 
fore, put 


«P( € ) = <P(€ . y,. P. T Q ) 


From the four quantities only one dimensionless com- 
bination can be formed: namely, 



The extent of the integration to 03 is merely a matter 
of form, $>( e ) decreases very substant ially with increasing 
e, and from a certain value of e on equals zero. 
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which, with z introduced as variable gives 


r , . 

i (a) lf,(0dz 


h = 


d 2g 


w 

J ° cp( z )ds 


/e\ K 

and with A = A !— ) gives 

0 A d / 


r 


h 


l y 2 
d 2g 


i 

J 


\ 0 
A o " 


Z n cp( z ) d z 


cp( z )d: 


y m pI n / 3 

jUi ni/'c 
u o 


= K 


1 V 2 u m D ffi/s 


d 2g d?‘r 0 m / 3 


(48) 


00 

o 

/ z cp( z )dz 


K = A 


0 co 


is a pure number. 


J 


cp( z )dx 


Considering the relation 


d rr 'Yhd 

"Yh drrT 0 i or T o = 

4 4 i 


existing between h and T 0 there is obtained 


h 


/ , ‘Y m / s d m / g = K J__ v® ® 

V s) 4 m/ £ l m 71 d W m 2g P 
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which, solved for h leaves 


h 


a 



m 


l ^ 

d 2g 


am 

V dv/ 


(49) 


This law corresponds exactly to the Bias ius resistance 
law for smooth pipes provided 


or 


2 m 
2 + m 


1/4 


m = 



By introducing a specific formula for the weight function cp - say, 
after the type of the law of error — the relationship Between 
the constants of Blasius* law for smooth pipes and the con- 
stants of the law of roughness can he ascertained. The writer 
hopes to be able to return to the further development of 
these arguments. 


If, on the other hand, 
troduced in (45) the result 


the law \ \ Q 

» a / 7 

with A Q e = 



i is 


is in- 


h = 





or 



According t,o the analogy between pipes and channels the 
velocity in a channel with gradient J and hydraulic radius 
P would be 

, , , 0.5 _ 0,64 

y = constant J 1 P 

*Mr. Prandtl states that he has arrived at the same 
result by'an entirely different process, 
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Incidentally, according to S, Manning, the empirical 

formula v = constant J 0,5 P 0 * 66 , according to lorchheimer , 

the formula v = constant J 0 * 5 P 0 * 7 , and according to 

Hermanek, the formula v = constant J 0,5 P 0, ° gives a good 
representation of the test data in rough channels (reference 
8, p. 70), 


Translation hy J, Vanier, 
National Advisory Committee 
for Aeronautics. 
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